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LONG-TERM  GOALS 

Key  to  assessing  the  risk  of  naval  activities  (such  as  sound  exposure)  on  marine  animals  is  an 
understanding  of  where  animals  occur  and  what  factors  motivate  these  movements.  The  rapid 
advancement  of  electronic  tracking  and  remote  sensing  technologies  has  enabled  researchers  to  link 
pelagic  predator  movements  and  oceanic  processes.  This  information  is  critical  for  understanding 
distribution  and  residence  time  of  vertebrates  within  an  ocean  area  and  for  managing  interactions  with 
anthropogenic  activities.  Marine  predators  interact  with  a  dynamic  ocean  that  change  on  time  scales 
ranging  from  minutes  to  millennia.  Knowledge  of  these  movement  interactions  is  incomplete  but 
critical  to  understanding  dynamic  distributions,  managing  anthropogenic  disturbance,  and  predicting 
responses  to  climate  change.  This  project  utilizes  the  largest  database  of  existing  marine  vertebrate 
tracking  and  behavior  data  to  build  upon  the  significant  advances  in  tag  technology,  data  analyses  and 
management  accomplished  under  the  Tagging  of  Pacific  Predators  (TOPP)  program.  This  is 
accomplished  by  establishing  a  behavioral  baseline  to  assess  the  potential  costs  of  displacement  in 
terms  of  reduced  foraging  success.  The  project  also  involves  a  synthesis  of  electronic  tracking  and 
remote  sensing  data,  focusing  on  a  cross-taxa  examination  of  marine  predator  distribution  and 
movement  patterns  to  identify  hotspots,  foraging  patterns  and  movement  corridors  in  the  California 
Current. 

OBJECTIVES 

1)  Identify  and  map  focal  feeding,  breeding,  and  migration  routes. 

2)  Model  spatio-temporal  oceanographic  habitat  utilization  and  predict  regions  of  animal 
occupancy  and  use  based  on  oceanographic  features. 
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3)  Utilize  this  model  framework  to  assess  the  impaet  of  displaeement  from  primary  feeding  areas 
due  to  disturbances  (such  as  acoustic  disturbance). 

APPROACH 

There  are  two  distinct  components  to  this  effort,  each  of  which  addresses  a  project  objective.  First,  we 
will  use  existing  TOPP  tracking  data  to  generate  overall  utilization  distributions  as  well  as  single 
species  distributions  and  further  categorize  track  segments  by  behavioral  state  using  a  combination  of 
state  spaced  models  and  the  fractal  landscape  method  to  determine  regions  of  area  restricted  search 
(ARS)  (Jonsen  et  al.  2003;  Jonsen  et  al.  2006;  Tremblay  et  al.  2007).  Next,  we  will  model  the  links 
between  oceanographic  parameters  and  animal  movement  patterns.  The  output  from  these  models  will 
be  used  to  develop  predictive  models  of  marine  vertebrate  distribution  based  on  oceanographic 
parameters. 

The  TOPP  data  management  system,  central  to  this  approach,  is  comprised  of  a  postgres  database, 
geospatially  enabled,  with  an  automated  tag  data  acquisition  and  processing  front  end,  using  a 
combination  of  in-house  produced  data  processing  tools  and  vendor  provided  tag  decoding  routines. 
Since  we  are  mainly  a  linux-based,  open  source  system,  and  the  vendor  tools  are  exclusively 
Windows-based,  the  vendor  tools  are  generally  run  in  a  virtual  Windows  machine  environment,  (using 
the  VirtualBox  software  package),  and  their  operation  is  automated  using  the  scripting  tool  Auto  It. 

Our  data  flow  starts  with  tag  data  input,  automatic  for  satellite  tag  data  and  semi-automatic  for  archival 
tag  data,  where  the  tags  must  be  manually  downloaded.  Once  tag  data  are  input  to  the  system,  they  are 
passed  through  our  QC  processing  system,  geolocated  where  appropriate,  and  loaded  into  the  database. 
Once  the  tag  information  is  in  the  database  we  can  provide  track  and  profile  data  in  combination  with 
oceanographic  parameters,  such  as  SST,  Chl-A,  and  sea  surface  height.  For  transmitting  tags,  this  can 
all  happen  in  near  real-time. 

WORK  COMPLETED 

Deployments  of  4,306  electronic  tags  yielded  1,791  individual  animal  tracks  from  23  species,  totaling 
265,386  animal  tracking  days  (Figure  1).  Technologies  had  different  location  precisions,  including: 
Argos  satellite  tags  (n=l,183),  archiving  and  satellite  transmitting  tags  (n=l,008),  and  archival 
geolocation  tags  (n=2,115)  the  latter  two  of  which  provide  estimates  of  position  based  on  sunlight, 
time,  and  sea  surface  temperatures  (SST).  These  different  positions  first  required  the  estimates  of  the 
precision  of  each  method,  coupled  with  incorporation  of  all  the  data  into  a  state  based  model.  Once 
this  was  accomplished  all  tracking  data  were  reduced  to  a  single  position  each  day  for  each  individual 
animal.  These  data  were  then  used  in  the  final  analysis. 
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Figure  1  \  All  TOPP  species  state  space  position  estimates  and  distribution  from  electronic  tagging, 
a,  Daily  mean  position  estimates  (circles)  and  annual  median  deployment  locations  (white  squares) 
of  all  tagged  species,  b,  Daily  mean  position  estimates  of  the  major  TOPP  guilds  (from  left):  tunas 
(yellowfin,  bluefin  and  albacore),  pinnipeds  (northern  elephant  seals,  California  sea  lions  and 
northern  fur  seals),  sharks  (salmon,  white,  blue,  common  thresher  and  mako),  seabirds  (Laysan  and 
black-footed  albatrosses  and  sooty  shearwaters),  sea  turtles  (leatherback  and  loggerhead)  and 
cetaceans  (blue,  fin,  sperm  and  humpback  whales). 


RESULTS 

The  tracking  data  reveal  that  the  California  Current  Large  Marine  Ecosystem  (CCLME; 
Supplementary  Eigure  1)  is  an  important  habitat  (Eigures  2-4)  for  tunas  (Pacific  bluefin,  Thunnus 
orientalis;  yellowfin,  T.  albacares),  sharks  (shortfin  mako,  Isurus  oxyrinchus;  white,  Charcharodon 
carcharias;  salmon,  Lamna  ditropis;  blue,  Prionace  glaucus',  common  thresher,  Alopias  vulpinus), 
cetaceans  (blue  whales,  Balaenoptera  musculus),  pinnipeds  (Northern  elephant  seals,  Mirounga 
angustirostris',  California  sea  lions,  Zalophus  californianus),  seabirds  (Eaysan  albatross,  Phoebastria 
immutabilis;  black-footed  albatross,  P.  nigripes;  sooty  shearwaters,  Puffinus  griseus),  and  turtles 
(leatherback  turtle,  Dermochelys  coriacea;  loggerhead  turtle,  Caretta  caretta).  Annual  migratory 
periodicity  was  evident  in  the  movements  of  many  tagged  animals  that  showed  fidelity  to  the  cool, 
nutrient  rich  waters  of  the  CCLME  (Eigure  2-3).  Extended  residency  within  the  CCLME  was  revealed 
by  examining  tracks  that  spanned  multiple  seasons  using  a  switching  state-space  model  (Eigure  3). 
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Numerous  species  exhibited  a  strong  attraction  to  the  CCLME  and  undertook  long  migrations  (>  2,000 
km)  from  the  western,  central,  or  south  Pacific  basin  (leatherback  sea  turtles,  black-footed  albatrosses, 
sooty  shearwaters,  bluefm  tuna,  and  salmon  sharks;  Figure  2b).  Species  exhibited  a  seasonally 
recurring  north-south  migration  (bluefin  and  yellowfin  tunas,  mako,  white,  and  salmon  sharks,  blue 
whales,  male  elephant  seals  and  leatherback  sea  turtles;  Figure  3g-l)  in  the  North  Pacific  and  within  the 
CCFMF.  Other  taxa  undertook  movements  between  nearshore  and  offshore  waters,  with  a  residency 
period  within  the  CCFMF  or  Gulf  of  Alaska,  followed  by  an  offshore  migration  that  ranged  into  the 
North  Pacific  Transition  Zone  (female  elephant  seals,  salmon  sharks,  Faysan  albatrosses),  the 
Subtropical  Convergence  Zone  (blue  and  mako  sharks,  leatherback  turtles),  or  the  eastern  Pacific  and 
Hawaiian  Islands  (white  sharks,  albacore  tuna,  and  black-footed  albatrosses).  The  mechanisms  and 
cues  underlying  fidelity  to  seasonally-modulated  migration  pathways  are  not  entirely  known,  but  may 
represent  a  capacity  to  discriminate  among  areas  of  seasonal  significance  for  foraging,  or  reproduction. 


Supplementary  Figure  1:  Schematic  of  dominant  ocean  current  features  in  the  north  Pacific  Ocean. 

The  California  Current  Large  Marine  Ecosystem  is  outlined  (dashed  line).  The  north  Pacific 

Transition  Zone  is  delineated  by  the  dotted  lines. 

The  CCFMF  is  a  highly  retentive  area  for  many  species  tagged  there,  and  an  attractive  area  for  animals 
undergoing  long  migrations  from  the  western  and  central  North  Pacific  and  Gulf  of  Alaska  (Figures  2; 
4a,b).  Pacific  bluefin  and  yellowfin  tunas  spent  significantly  more  time  in  the  CCFMF  than  expected 
based  on  null  model  simulations.  Several  species  -  tunas;  lamnid  sharks;  leatherback  turtles;  and  blue 
whales  had  substantial  residency  periods  within,  or  were  return  migrants  to  the  CCFMF  based  on 
behavioral  inferences  from  a  switching  state-space  model  (Figure  4c).  Additionally,  all  species  tagged 
outside  the  CCFMF  spent  significantly  more  time  on  average  in  the  CCFMF  than  expected  based  on 
null  model  simulations.  Representatives  from  several  guilds  exhibited  cross-basin  migrations  (>  2,000 
km)  into  the  CCFMF  from  the  western  (leatherback  turtle,  bluefin  tuna),  central  (black-footed 
albatross,  salmon  shark),  and  south  Pacific  basins  (sooty  shearwaters;  Figure  2b). 
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Figure  2  \  Fidelity  and  attraction  to  the  CCLME.  Left,  Examples  of  pelagic  predators  released  and 
electronically  tracked  in  the  CCLME  that  show  fidelity  to  deployment  locations  and  the  CCLME. 

We  show  the  release  locations  (square),  pop-up  satellite  end  point  locations  (triangle)  and  daily 
mean  positions  (circles)  of  the  following  species:  yellowfin  tuna  (yellow),  bluefin  tuna  (white),  white 
shark  (red),  elephant  seal  (blue)  and  salmon  shark  (orange).  Right,  Individual  tracks  of  pelagic 
animals  released  .2,000km  away  from  the  CCLME  that  are  indicative  of  cross-basin  or  ecosystem 
attraction  to,  and  temporary  residency  within,  the  eastern  North  Pacific.  Symbols  are  as  in  a,  for 
leatherback  sea  turtles  (green),  sooty  shearwaters  (pink),  fur  seals  (pale  yellow),  black-footed 

albatrosses  (black)  and  salmon  sharks  (orange). 
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Eigure  3  \  Latitudinal  migration  cycles  and  seasonal  climatologies  within  the  CCLME.  a.  Monthly 
mean  latitudes  of  predators  residing  within  or  migrating  to  the  CCLME.  Black  line  segments  denote 
gaps  where  no  data  were  available.  Sample  sizes  indicate  the  numbers  of  individual  tracks 
contributing  to  the  time  series,  b.  Seasonal  climatologies  in  the  California  Current  for  tunas  (Pacific 
bluefin,  blue;  yellowfin,  black),  sharks  (salmon  shark,  brown;  shortfin  mako,  black;  white  shark, 
blue)  and  blue  whales  relative  to  median  chlorophyll  a  densities  and  SST  values 

between  2000  and  2009. 
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The  retention  with  and  attraction  to  the  CCLME  is  consistent  with  the  high  productivity  of  this  region 
that  supports  large  biomasses  of  krill,  sardines,  anchovies,  salmon,  groundfish,  and  squid  providing  a 
predictable  forage  base  for  top  predators.  The  NPTZ  is  another  important  region,  serving  as  an  east- 
west  migration  corridor  (Figure  4a)  and  foraging  region  for  elephant  seals,  salmon  and  blue  sharks 
(Figure  4c),  Faysan  and  black-footed  albatrosses,  and  bluefm  tuna  (Figure  1).  This  is  a  complex  region 
encompassing  an  abrupt  north-to-south  transition  between  subarctic  and  subtropical  water  masses  with 
dynamic  frontal  regions. 


Figure  4  \  Predator  density  maps  and  residency  patterns,  a,  Density  of  large  marine  predators  within 
the  eastern  North  Pacific.  Densities  of  the  time-weighted  and  species-normalized  position  estimates 
of  all  tagged  individuals  were  summed  within  lu31u  grid  cells,  b,  Density  of  large  marine  predators 
within  the  CCLME  at  a  0.25u30.25u  resolution,  c,  Patterns  of  resident  (slow,  arearestricted 
movements)  versus  transient  (fast,  directed  movements)  behaviours  of predators  that  primarily 
occupied  or  migrated  to  the  CCLME,  estimated  using  a  switching  state-space  model.  The  coloured 
points  grading  fromblue  to  yellow  display  the  posterior  mean  probability  of  the  resident  behaviour 
associated  with  each  daily  mean  position  estimate.  Each  panel  displays  residency  patterns  for  ten 
individuals.  Uncertainty  in  position  estimates  in  a  and  b  is  included  by  calculating  densities  using  all 
2000  Markov  chain  Monte  Carlo  samples  from  the  joint  posterior  distributions  of  the  daily  positions, 
rather  than  using  only  the  posterior  means.  SST  contours  in  a  are  denoted  by  solid  white  lines. 

Exclusive  economic  zones  are  delineated  by  solid  black  lines. 

The  second  component  of  this  effort  was  to  identify  the  oceanographic  parameters  responsible  for  the 
distribution  of  these  predators.  To  investigate  which  aspects  of  the  biophysical  environment  putatively 
attract  these  predators,  we  explored  both  presence/absence  and  relative  habitat  use  with  generalized 
additive  mixed  models.  We  examined  the  collective  response  of  16  marine  predator  species  to 
environmental  covariates.  In  the  binary  presence/absence  model,  predator  incidence  showed  a  strong 
positive  relationship  with  SST,  across  a  broad  temperature  range  and  peaking  near  15°C.  Tagged 
animals  occupied  a  small  portion  of  cool,  nutrient-rich  water  in  coastal  regions  and  northern  latitudes 
compared  to  the  broadly  available  warm  oligotrophic  waters  in  lower  latitudes.  The  strong  positive 
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relationship  between  relative  density  and  Chl-a  suggests  the  suite  of  tagged  speeies  are  preferentially 
occupying  regions  of  high  productivity.  The  observed  patterns  of  predator  distribution  in  this  study 
may  be  indicative  by  trade-offs  between  thermal  tolerances,  either  directly  by  the  predators  or 
indirectly  on  lower  trophic  levels,  and  access  to  areas  of  higher  productivity. 

To  examine  how  closely  related  taxa  partition  marine  resources,  we  compared  thermal  preferences 
from  in-situ,  tag-based  SST  measurements  for  sympatrically  occurring  species  within  three  guilds 
(albatrosses,  tunas,  sharks).  Differences  in  habitat  utilization  evident  among  congeneric  species 
illustrate  how  recently  divergent  species  partition  the  oceanic  environment  (Figure  5).  During  the  June 
to  November  post-breeding  phase,  black-footed  albatrosses  were  associated  with  a  broader  and  warmer 
range  of  SST  primarily  in  the  eastern  Pacific,  whereas  Laysan  albatrosses  were  associated  with  a 
narrower  and  colder  range  of  SST  in  the  western  and  central  North  Pacific  (Figure  5a).  Bluefin  tuna 
ranged  farther  north  in  the  colder  waters  of  the  CCLME,  whereas  yellowfin  occupied  warmer  waters  in 
the  southern  CCLME  (Figure  5b).  These  differences  are  consistent  with  physiological  specializations 
in  bluefin  tuna  cardiac  performance.  The  lamnid  sharks  had  a  more  complex  separation  of  habitats. 
Salmon  sharks,  with  their  cold  tolerant  cardiac  physiology,  occupied  the  cooler,  subarctic  waters  in  the 
North  Pacific.  Most  salmon  sharks,  but  not  all,  migrated  seasonally  into  the  warmer  NPTZ  and 
CCLME  waters  (Figure  5c).  White  sharks  overlapped  with  salmon  sharks  in  the  nearshore  CCLME, 
but  also  migrated  into  warmer,  offshore  waters  of  the  Subtropical  Gyre  and  the  Flawaiian  Islands 
(Figure  5  c).  Shortfin  makos  were  distributed  throughout  the  CCLME  and  into  the  Subtropical  Gyre, 
but  occupied  a  thermal  range  intermediate  to  the  two  modes  of  the  white  shark  range  (Figure  5  c). 


Figure  5  \  Niche  separation  within  three  predator  guilds.  Spatial  distribution  and  thermal 
habitat  use  (insets)  across  three  guilds  of  sympatric  species:  Laysan  and  black-footed 
albatrosses  (n555  individuals,  8,743  daily  SSTs;  a).  Pacific  bluefin  and  yellowfin  tunas 
(n5376  individuals,  75,1 77  daily  SSTs;  b),  lamnid  (salmon,  mako  and  white)  sharks 
(n5137  individuals,  12,971  daily  SSTs;  c).  SST  profiles  are  daily  means  of  tag-derived  SSTs. 
Linear  mixed-effects  model  estimates  of  mean  (695%  confidence  interval)  SST  for  each 
species  are  displayed  as  diamonds  at  the  top  of  each  inset  graph. 
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Figure  6.  An  energy  landscape  for  Pacific  Bluefin  Tuna  between  March  and  May  2004. 

Warm  colors  represent  areas  with  higher  energy  acquisition,  while  cooler  colors  represent 
areas  with  lower  energy  acquisition  (kcals/day).  The  overlaid  black  dots  show  daily  positions. 

To  study  the  ability  to  use  eleetronie  tagging  to  investigate  building  energetie  landseapes  in  the  wild 
we  ran  validation  experiments  in  a  eaptive  feeding  study.  Bluefin  tuna  warm  their  viseera  during 
digestion.  This  warming  is  thought  to  be  due  to  a  eombination  of  faetors,  foremost  among  whieh  is 
speeific  dynamie  aetion  (SDA),  the  heat  produeed  by  digestion,  absorption  and  anabolie  proeesses 
within  the  stomaeh  and  eaeeal  mass.  The  duration  of  the  elevated  viseeral  temperature  is  strongly 
eorrelated  with  the  ealorie  value  of  an  ingested  meal.  We  have  eondueted  an  experiment  using 
arehival-tagged  eaptive  bluefin  tuna  in  the  Tuna  Researeh  and  Conservation  Center  to  quantify  the 
magnitude  of  viseeral  warming  assoeiated  with  differing  ration  sizes  of  squid  and  sardine,  at  different 
ambient  temperature  levels.  The  data  eolleeted  from  these  experiments  were  used  to  develop  a  model 
whieh  ean  be  applied  to  arehival  tag  data  eolleeted  from  wild  Paeifie  bluefin  tuna,  in  order  to  estimate 
daily  energy  intake.  This  measure  of  foraging  sueeess  is  now  being  used  to  identify 
environmental/oeeanographie  eorrelates  of  foraging  sueeess  as  well  as  eharaeterizing  foraging 
“hotpots”  for  Paeifie  bluefin  tuna.  So  far.  Generalized  Additive  Models  (GAMs)  have  been  used  to 
produee  seasonal  energy  maps  for  the  California  Current,  highlighting  areas  where  daily  energy  intake 
is  highest  (Figure  6). 

IMPACT/APPLICATIONS 

Critieal  to  determining  the  impaet  of  exposure  to  naval  operations  on  marine  animals  is  relating  the 
intensity  and  duration  of  an  exposure  to  the  time  animals  spend  in  proximity  to  the  souree,  and  the 
biologieal  funetion  of  that  time.  The  proposed  predietive  models  of  eritieal  marine  animal  habitat 
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utilization  are  the  essential  behavioral  components  to  determine  whether  and  where  naval  operations 
might  impact  marine  mammals  and  other  marine  vertebrates. 
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